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?H ■ ABSTRACT 

> 

. Context. The Red MSX Source (RMS) survey has identified a large sample of candidate massive young stellar objects (MYSOs) and ultra 
compact (UC) HII regions from a sample of ~2000 MSX and 2MASS colour selected sources, 
j Aims. To search for HtO masers towards a large sample of young high mass stars and to investigate the statistical correlation of H2O masers 
^ with the earliest stages of massive star formation. 

Methods. We have used the Mopra Radio telescope to make position-switched observations towards ~500 UCHII regions and MYSOs candi- 
dates identified from the RMS survey and located between 190° < / < 30°. These observations have a 4<r sensitivity of ~1 Jy and a velocity 
£" — . . resolution of ~0.4 km s _1 . 

r- H Results. We have detected 163 H2O masers, approximately 75% of which were previously unknown. Comparing the maser velocities with the 
velocities of the RMS sources, determined from 13 CO observations, we have identified 135 RMS-H 2 maser associations, which corresponds 
to a detection rate of ~27%. Taking into account the differences in sensitivity and source selection we find our detection rate is in general 
agreement with previously reported surveys. 
f^**) , Conclusions. We find similar detection rates for UCHII regions and MYSOs candidates, suggesting that the conditions needed for maser activ- 
L" ity are equally likely in these two stages of the star formation process. Looking at the detection rate as a function of distance from the Galactic 
. centre we find it significantly enhanced within the solar circle, peaking at ~37% between 6-7 kpc, which is consistent with previous surveys of 
^ . UC HII regions, possibly indicating the presence of a high proportion of more luminous YSOs and HII regions. 
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1. Introduction 

Massive young stellar objects (hereafter MYSOs) are in the 
early phase in the life of OB stars when fusion has most likely 
started in the core, but when they have not yet begun to ion- 
ize their surroundings to form an HII region. The embedded 
mid-infrared poi nt source usually possesses a strong ionized 



stellar wind (e.g.. [Bunn et al 
lar outflows (e.g.. lLadalll985 



1995) and drives bipolar molecu- 



Bronfman et al ] |2008t and refer 



ence therein). This brief (10 4 ~ 5 years) phase is clearly crucial 
to our understanding of how these massive young stars form 
since during this time any ongoing accretion will be halted and 
the final mass of the star is set. In addition, the winds, outflows 
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and eventual HII regions have an important feedback role in 
determining the fate of the rest of the molecular cloud. 

The main difficulty in studying MYSOs is their relative rar- 
ity. Most current work still relies on samples drawn from the 
IRAS Point Source Catalogue, but whose selection criteria bias 
them away from complex regions by requiring them to be away 
from HII regions identified in single-dish r adio surveys (e.g., 
Molinari et al]|l996t ISridharan etail 12002). Additionally, the 



known samples are too small to test many aspects of massive 
star formation theories, and are probably unrepresentative of 
the class as a whole. In an effort to address these issues we 
have conducted a new galaxy-wide search for MYSOs start- 
ing from the MSX mid-infrared survey of the Galactic Plane 



([Price et al.ll2001l) . This has significantly better spatial reso- 



lution than IRAS and so does not have the same confusion 
problems. We ha ve used the MSX Point Source Catalogue 
(lEgan et al.ll2003b to colour-select a large sample of candidate 
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MYSOs ( Lumsden et al.ll2002 ). This initial sample was further 
refined using 2MASS data to eliminate blue objects and from 
visually inspecting the MSX images to remove move extended 
sources. Our colour-selection and subsequen t filtering pro- 
duced a sample of ~2000 MYSO candidates dLumsden et al 
2002h. 



The Red MSX Source (RMS) Survey is a multi-wavelength 
programme of follow-up observations designed to identify 
genuine MYSOs and ultra compact (UC) HII regions from 
the other kinds of embedded or dusty objects such as plan- 
etary nebulae (PNe), evolved stars and nearby low-ma ss 



YSOs (IHoare et alj|2005b iMottram et all 120061: lUrquhart et al 
2008bl) . The main aim of the project is to produce a large 
sample of MYSOs and UCHII regions, and to compile a 
database of complementary multi-wavelength data with which 
to study their properties^] We have used arcsecond resolu- 
tion mid-infrared ima ging from the Spitzer GLIMPSE survey 
(Benjamin et al ]|2003l) or our own ground-based imaging (e.g., 
Mottram et al .120071) to reveal multiple and/or extended sources 
within the MSX beam, as well as MYSOs in close proximity 
to existing HII regions. We have obtained arcse cond resolution 



radio continu um data with ATCA and the VLA (lUrquhart et al 



2007al 120091) to identify UCHI I regions and PNe, whilst oh 
servations of 13 CO transitions (lUrquhart et al.l l2007bl l2008al) 
deliver kinematic velociti es. Finally we hav e obtained near- 
infrared spectroscopy (e.g. JClarke et al hood) which allows us 
to identify the more pathological evolved stars. 

The kinematic velocities can be used in conj unction 



with a Galactic rotation model (e .g., iBrand & Blita 11993 
Alvarez et al. Il990b Clemens 1985) to derive kinematic dis 
tances and luminosities, which allow us to distinguish be- 
tween nearby low- and intermediate-mass YSOs and genuine 
MYSOs. However, the velocity of sources located within the 
solar circle - which accounts for ~80% of the sample - results 
in two possible kinematic distances equally spaced on either 
side of the tangent position; these are referred to as the near 
and far distances. This distance ambiguity needs to be solved 
before luminosities can be calculated and our sample of MYSO 
candidates can be turned into a sample of bona fide MYSOs. 
We have solved the distance ambiguities towards a sample of 
MYSOs using the HI self-absorption (see iBusfield et alil2006l 
for details and a description of the technique) but distances to 
the whole sample are not currently available. 

By combining these observational data sets it is possible to 
identify the different source populations and remove the con- 
taminates. With the classification effectively complete we are 
now moving into the exploitation phase of the RMS survey. 
The first step is to examine the global characteristics of this 
galaxy-wide sample of massive young stars as a prelude to de- 
tailed studies of sub-samples. Part of this involves determining 
the physical and chemical nature of the environment as a way 
of gauging their evolutionary status. We have used the Mopra 
telescope and the large bandwidth available with the broad- 
band (~8 GHz) MOPS spectrometer to survey a sub-sample of 
MYSO candidates and UCHII regions at 12 mm. The primary 
focus is to study water maser emission and hot ammonia, but 

1 http://www.ast.leeds.ac.uk/RMS 



also class II methanol maser emission and other serendipitous 
molecular transitions found in the range. 

In this paper we present the results of our search for water 
masers made towards ~500 MYSO candidates and UCHII re- 
gions observable by Mopra. Water masers are the most likely 
masing species to be found associated with the mid-infrared 
bright MYSO phase. VLBI proper motion measurements show 
that the y arise in the main from the base of bipo lar outflows 
and jets dGoddi et alj2005tlMoscadelli et al.l2005l) and as such 
trace an important part of the physics of massive star forma- 
tion. Although the spatial resolution is not sufficient to warrant 
a detailed study of individual sources, these single-dish obser- 
vations allow us to look at the statistical occurrence of water 
masers towards massive star forming regions, as well as paving 
the way for future interferometric follow-up observations. 



2. Observations and source selection 

2.1. Source selection 

To date our ongoing programme of classification has led to the 
identification of approximately 500 MYSO candidates and a 
further 600 UC HII regions. In addition to these categories we 
have two further classifications that are relevant to this study; 
these are the HII/YSO and Young/old classifications. The first 
of these is reserved for RMS sources where both an HII region 
and YSOs are found within the MSX beam (i.e., <18")- We 
have identified ~200 of these, which are interesting in that they 
may indicate sequential star formation. The second of these 
classes contains sources that appear isolated in near-infrared 
images and which would normally indicate an evolved star 
classification, however, they also display strong CO towards 
them and/or lie on dark filaments seen at 8 /an. They could be 
luminous evolved background stars seen through a molecular 
cloud or be genuine YSOs forming in the cloud. These cannot 
be unambiguously classified until infrared spectra have been 
obtained and so remain within this class until a reliable classi- 
fication can be attributed. 

Approximately 900 RMS sources classified as one of the 
four types mentioned in the previous paragraph are observable 
by the Mopra telescope (i.e., 190° < I < 30°). However, due to 
a combination of limited observing time and density of sources 
towards the Galactic centre we were only able to observe a lit- 
tle over half this number. Since the overall aim of the RMS 
project is to identify a large sample of massive YSOs, prior- 
ity was given to sources identified as either YSOs or HII/YSOs 
with ~80% and ~9Q% of these types observed respectively. A 
high priority was also given to the Young/old class of sources, 
with ~80% being observed, as a significant number of gen- 
uine YSOs may still be included in this catagory. Many of the 
sources excluded from our observations of these three classes 
of objects were on the grounds that water masers have previ- 
ously been reported in the literature, or because they are located 
within a region of the Galactic Plane which is already the fo- 
cus of a blind 19- 27 GHz survey currently underway (HOPS; 
Walsh et al.ll2008l) . HII regions were considered a lower prior- 
ity and were observed when time and RA range became avail- 
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Table 1. Summary of sources observed, positions and sensitiv- 
ities. 



MSX name 


RA 


Dec 


r.m.s 




J2000 


J2000 


(mK) 


G018. 8330-00.3004 


18:26:23.64 


-12:39:40.3 


23 


GO 19.7540-00. 1279 


18:27:31.48 


-11:45:56.8 


25 


G027.9334+00.2056 


18:41:34.15 


-04:21:11.8 


14 


G289. 1447-00.3454 


10:57:07.24 


-60:07:15.5 


26 


G293.5607-00.6703 


11:30:07.10 


-62:03:12.9 


24 


G296.8926-01.3050 


11:56:49.89 


-63:32:04.9 


24 


G300.7221+01.2007 


12:32:50.83 


-61:35:28.6 


23 


G301. 8 147+00.7808 


12:41:53.78 


-62:04:13.0 


26 


G305. 9402-00. 1634 


13:17:52.94 


-62:52:50.5 


25 


G311.5131-00.4532 


14:05:45.91 


-62:04:49.4 


25 


G3 11.5671+00.3189 


14:04:22.27 


-61:19:26.7 


26 


G3 11.9799-00.9527 


14:10:51.43 


-62:25:16.3 


20 


G313.5769+00.3267 


14:20:08.35 


-60:41:55.3 


23 


G319.8366-00.1963 


15:06:54.88 


-58:32:57.4 


23 


G320.2878-00.3069 


15:10:19.60 


-58:25:07.3 


25 


G320.3767-0 1.9727 


15:17:39.02 


-59:47:49.2 


23 


G325.3566-00.0422 


15:40:11.95 


-55:23:19.3 


20 


G327.9205+00.0921 


15:53:32.92 


-53:42:05.4 


21 


G328.2658+00.5316 


15:53:27.93 


-53:08:33.7 


21 


G333.3752-00.2015 


16:21:05.95 


-50:15:14.0 


21 


Notes: Only a small 


portion of 


the data is provided 



the full table is only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A+ A/ 



able; only approximately 30% of the available HII regions have 
been observed. 



2.2. Observations and data reduction 

The observations were made using the Mopra 22m 
Radiotelescope in April and September 2008 towards 
499 massive star forming regions. Mopra is located near 
Coonabarabran, New South Wales, Australia^ The telescope 
is situated at an elevation of 866 metres above sea level, and at 
a latitude of 3 1 degrees south. 

The telescope is equipped with a 12 mm receiver with a fre- 
quency range of 16 to 27.5 GHz. The UNSW Mopra spectrom- 
eter (MOPS) is made up of four 2.2 GHz bands which over- 
lap slightly to provide a total of 8 GHz continuous bandwidth. 
Up to four zoom windows can be placed within each 2.2 GHz 
band allowing up to 16 spectral lines to be observed simultane- 
ously. Each zoom window provides a bandwidth of 137 MHz 
with 4096 channels, which at 12 mm corresponds to a total 
velocity range of ~2000 km s with a velocity resolution of 
~0.4 km s _1 per channel. 

The 8 GHz bandpass was centred at 23.5 GHz provid- 
ing complete coverage over the 19.5-27.5 GHz range with a 
corresponding resolution of 2.7-1.8' respectively. Individual 
zoom windows were deployed to cover: the water maser at 
22.235 GHz, the NH 3 (1, 1) and (2, 2) inversion lines at 

2 Mopra is operated by the Australia Telescope National Facility, 
CSIRO. 
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Fig. 1. Histogram of the r.m.s. noise level in the different ob- 
servations. The bin size is 5 mK channel -1 . 



23.694 and 23.722 GHz, two more covering the NH 3 (3, 3) 
and (4, 4) lines located at 23.870 and 24. 139 GHz respectively. 
Additional zoom windows were placed to cover the thermal 
methanol lines around 24.9 GHz, and at 23.121 GHz to cover 
the rare methanol maser transition (Cragg et al. 2004). The re- 
maining windows were used to search for the more common 
19.97 GHz methanol maser line (Ellingsen et al. 2004) and 
HC5N (7-6) transition at 18.64 GHz. In this paper we will con- 
centrate on the detection statistics of water masers and post- 
pone the analysis of the other transitions to a subsequent pub- 
lication. 

The observations were conducted in position-switched 
mode towards 499 RMS sources (see Table [T] for positions). 
Each source was observed for a total of ten minutes of on- 
source integration, split into a number of separate scans con- 
sisting of 1 minute on- and 1 minute off-source. Reference po- 
sitions were offset from the MSX positions by 1 degree in a di- 
rection perpendicular to the Galactic plane. Given that the size 
of the beam at 12 mm (~2') is much larger than the possible 
deviations from the global pointing model no pointing checks 
were performed. System temperatures were between 65-100 K 
depending on weather conditions and telescope elevation, re- 
sulting in typical rms values of ~25 mK per channel (see Fig.Q] 
for distribution). 

Cross scans of Jupiter were made to determine the telescope 
main beam efficiency during our observations which was found 
to be 77MB ~0.55 at 23 GHz. The corrected antenna tempera- 
tures (T* A ) were put on the main beam temperature scale by di- 
viding by the telescope main beam efficiency (77mb)- Finally, the 
main beam temperatures were converted to flux density using 
a conversion factor of 6.41 Jy K 1 (Eqn 7.19: lRohlfs & Wilson! 
20041) . Taking the conversion from T^ to Jy into account the 
nominal Act detection limit is ~1 Jy per channel. Hereafter all 
quoted intensities will be given in Janskys. 

The data were reduced using the ATNF spectral line reduc- 
ing software (ASAP). The reduction steps consisted of divid- 
ing the individual on-off scans to remove sky emission (these 
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Fig. 2. Spectra of detected HhO masers. Only a small portion of the plots are provided here, the full fig- 
ure is only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat 7J/A+A/ 



scans were then inspected and poor scans were removed), fit- 
ting a low-order polynomial to the baseline. Finally, the scans 
were averaged together to produce a single spectrum for each 
source. 



3. Results and analysis 

3.1. Detection statistics and general properties 

We detect maser emission towards 153 of the 499 star form- 
ing regions observed. In 10 of these cases two distinct groups 
of maser emission were observed separated in velocity by > 
30 km s _1 . The emission structure and velocity separation 
makes it likely that the emission seen in these cases is from two 
distinct sites of maser emission along the same line of sight, 
within the ~2 arcmin beam. We have therefore classified these 
as separate detections and indicate these in the plot presented 
in Fig. [2] and Table [2] by appending an to the MSX name 
(where n = 1 or 2). In all of these cases at least one of the line 
of sight components is found at a similar velocity as a molecu- 
lar tracer and is therefore likely to be associated with an RMS 



source (see Sect. 4.1 for details). However, the origin of these 
other masers detected along each line of sight is unclear since 
the majority (8) have velocities outside the velocity range usu- 
ally associated with molecular gas and so were not covered by 
our CO observations (see discussion in Sect. 4.3). Only two of 
these masers are found within the velocity range of our CO ob- 
servations. One of which is found to have a similar velocity as 
a weak CO component, whilst no CO emission is seen around 
the velocity of the other, however, given the difference in res- 
olution this is per haps not surprising (C O observations have a 
resolution of 34"; IUrquhart et al.|[2007bb . 

Taking these cases into account, we have detected a total 
of 163 masers above a 4<x noise level (~1 Jy), which corre- 
sponds to an overall detection rate of ~32%. We conducted a 
literature search using SIMBAD and found that approximately 
75% of these masers were previously unknown. In Fig. [2] we 
present plots of all the detected masers. Where two masers have 
been detected along the same line of sight we present a separate 
plot for each. To avoid confusion the velocity range for a num- 
ber of these multiple detections has been reduced. As can be 
seen from the spectra presented in this figure, the maser emis- 
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Table 2. Parameters of detected masers. 



MSX Name 


RA 


Dec 


Peak Flux 


Vlsr 


V m in 


Vmax 




(J2000) 


(J2000) 


(Jy) 


(km s ') 


(km s ') 


(km s ') 


G010. 3844+02. 2128 


18:00:22.68 


-18:52:08.0 


8.4 


-5.8 


-9.1 


-2.5 


G010. 4413+00. 0101 


18:08:38.23 


-19:53:57.4 


2.8 


72.7 


63.5 


84.9 


G010. 5067+02. 2285 


18:00:34.58 


-18:45:17.6 


3.4 


20.8 


17.4 


22.1 


G011.4201-01.6815 


18:16:56.86 


-19:51:07.2 


4.4 


0.5 


-1.4 


1.6 


G012. 1993-00.0342 


18:12:23.66 


-18:22:51.6 


4.0 


2.1 


-0.2 


30.0 


G013. 0105-00. 1267 


18:14:22.25 


-17:42:47.8 


2.4 


5.3 


-7.4 


16.4 


G014. 9958-00.6732 


18:20:19.43 


-16:13:31.0 


18.9 


23.9 


3.6 


63.6 


G015.0939+00.1913 


18:17:20.86 


-15:43:47.2 


8.1 


30.2 


23.4 


33.8 


G015. 1288-00.67 17_1 


18:20:34.75 


-16:06:26.2 


1.7 


-41.9 


-43.0 


-41.3 


G015. 1288-00.6717.2 


18:20:34.75 


-16:06:26.2 


0.7 


22.2 


19.7 


25.6 


G016.8055+00.8149 


18:18:25.96 


-13:55:37.9 


5.1 


70.4 


68.0 


72.5 


GO 16.9270+00.9599 


18:18:08.59 


-13:45:05.7 


1 .7 


-65.8 


-68.0 


-59.1 


G017.9789+00.2335 


18:22:49.08 


-13:09:59.0 


4.7 


18.1 


15.8 


20.3 


G018. 1409-00.3021 


18:25:04.46 


-13:16:26.7 


0.7 


13.0 


10.7 


13.6 


G024.2 138-00.0439 


18:35:35.80 


-07:46:22.4 


4.0 


55.7 


53.2 


56.8 


G025. 4018+00. 0198 


18:37:34.19 


-06:41:18.6 


7.9 


-13.3 


-15.9 


-6.5 


G026. 1 094-00.0944 


18:39:17.14 


-06:06:43.9 


16.4 


25.1 


20.6 


28.6 


G028. 3271+00. 1617 


18:42:26.90 


-04:01:23.8 


1.9 


31.3 


29.1 


32.8 


/"■mn a 1 1 7 nn 7777 
OU 1 1 / — uu.zz / / 


15.4 / . JV.U4 


— UZ.ZU.4 /. / 


7 q 
Z.o 


1 o 


Q7 7 


1 n7 n 


G030.8667+00.1141 


18:47:15.81 


-01:47:08.8 


2.9 


37.8 


36.3 


48.4 


G188.8120+01.0686 


06:09:17.85 


+21:50:49.9 


10.4 


-5.9 


-8.3 


-3.8 


G200.0789-0 1.6323 


06:21:47.85 


+ 10:39:20.8 


2.0 


35.2 


33.8 


42.0 


G224.3494-02.0143 


07:05:12.43 


-11:04:31.0 


3.8 


16.4 


14.6 


21.4 


G224.6573-02.4978 


07:04:01.65 


-11:34:12.3 


0.6 


11.5 


9.8 


13.2 


G226.2728-00.4633 


07:14:26.80 


-12:03:58.3 


2.9 


17.3 


15.6 


17.9 


G229.571 1+00. 1525 


07:23:02.08 


-14:41:32.2 


9.0 


60.8 


47.7 


62.4 


G234.5093-00.3181 


07:31:08.76 


-19:15:29.8 


1.6 


43.5 


42.0 


45.4 


G236.8158+01.9821 


07:44:27.93 


-20:08:33.0 


15.2 


54.4 


39.4 


64.3 


G240.3156+00.0715 


07:44:52.00 


-24:07:41.1 


18.3 


65.3 


62.5 


78.9 



Notes: Only a small portion of the data is provided here, the full table is only available in electronic form at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.125.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



sion generally consists of either a single component or a group 
of emission peaks spread over a clearly identifiable velocity 
range. In order to include all emission components associated 
with each maser the plots presented in Fig. [2]cover a velocity 
range of ~100 km s _I and are centred on the velocity of the 
strongest component. 



In Table 12 we present the parameters of the detected H2O 
masers; we give the MSX name and J2000 co-ordinates in 
Columns 1-3. In Column 4 and 5 we give the peak flux density 
and the velocity of the peak emission. We give the minimum 
and maximum velocity range of the emission in Columns 6 
and 7. 



In Fig. [3] we present histograms of the distribution of the 
peak flux density and maser group velocity dispersion (defined 
as a number of discrete components within 30 km s of each 
other). The distribution of peak flux densities spans a range of 
4 orders of magnitude from ~0.5 Jy up to a couple of thousand 
Jy, however, more typical values are between 1-10 Jy. The his- 
togram of the distribution of maser velocity dispersion, shown 
in the right panel of Fig. [3] illustrates the large spread in ve- 
locity of maser emission. Maser groups consisting of a single 
component have an intrinsic width of as little as a few km s 
with multiple component groups having velocity spreads of up 
to 90 km s . However, these are the extremes with the major- 
ity of water masers having velocity ranges of 10-15 km s . 
This large velocity dispersion is rather typical for water masers 
and is one of the reas ons a connection with molecular o utflows 
is often inferred (e.g., Elitzur et al.lll989 lMentenlll996l) . 



4. Discussion 

4.1. RMS-maser associations 

We have detected and identified 163 water masers towards 
RMS sources. Not all will necessarily be associated, however, 
as some are likely to be chance alignments along the line of 
sight. Given the spatial resolution of these observations, we are 
unable to confirm associations but we can identify likely RMS- 
H2O associations by considering the velocities of the maser 
and the RMS sources (determ ined from CO observations, see 
Urquhart et al .l2007bl l l2008al for details) and look for a velocity 
correlation. 

In an effort to identify likely RMS-FbO maser associa- 
tions we compared the peak maser velocity with the velocity 
of the RMS source along each line of sight. The distribution 
shown in Fig. [4] illustrates the excellent correlation between 
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Fig. 3. Histograms of the H2O maser parameters. In the left panel we present the histogram of the peak intensity distribution (the 
bin size is 0.5 dex) and in the right panel we present a histogram of the H2O maser groups dispersions (binned using a value of 5 
km s ). 



the velocity of the peak intensity component in each maser 
group and the velocity of the molecular cloud within which 
the RMS source resides. The distribution is strongly peaked at 
< 2 km s , and falls off steadily until reaching a background 
level at ~ 14 km s . The shape of the histogram corresponds to 
that of a one-sided Gaussian and we have therefore determined 
the standard deviation from a Gaussian fit to the profile. 

We consider an RMS source and H2O maser to be associ- 
ated if the difference in their velocities is less than 3<x, where cr 
is the standard deviation of th e distribution ~5 km s" 1 (cf. a me- 
dian difference of 4.5 km s _1 ; Kurtz & Hofnerl2005 ). Applying 
this criterion (i.e., AV LSR < 15 km s -1 ) we find 131 RMS- 
H2O maser associations. No velocity has been assigned to five 
RMS sources due to the presence of multiple CO components 
in the sp ectra, which makes it d ifficult to determine a unique 
velocity ( Urquhart et al. 2007bh . However, it is still possible 
that the masers detected towards these sources are associated. 
To investigate these possible associations we compared the ve- 
locities of the peak maser components to the velocities of the 
detected CO components seen towards all five RMS sources. 
We found the peak component velocity of the maser within 
15 km s _I of at least two of the CO peaks present in four of 
the five cases, and therefore consider these to be possible as- 
sociations. Unfortunately, we were not able to use the presence 
of the masers to identify the CO component associated with the 
RMS source in these cases. Including the 131 RMS sources al- 
ready found to have an associated water maser this brings the 
total to 135. The remaining 27 masers are likely to be associ- 
ated with other star forming regions that happen to be located 
along the same line of sight as the RMS source. 



4.2. Comparison of RMS-H 2 detection rate 

Our detection rate for finding water masers associated with 
RMS sources is ~27 + 2%. Comparing this value with other 
detection rates reported from previous H2O surveys we find 
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Fig. 4. Histogram of the differences between the velocity of the 
detected H2O maser and the 13 CO velocity of the RMS source 
within the beam. The bin size is 2 km s -1 . 



it is towards the lower end of a r ange of values (i. e ., 20- 
67 %). For example, observ ations by Churchwell et al. ( 1990l) 
and iKurtz & Hofnerl (12005b towards bright IRAS sources re- 
ported relativ ely high de t ection rates of 67% and 5 5 % res pec- 
tively, while IPalla et all dl99ll) and ICodella et all d 19951) re- 
ported detection rates of only 26% and 20% respectively. 

All of th ese surveys where made tow ards IRAS sources that 
satisfied the Wood & Churchwell J 19891) colour selection crite- 
ria for identification of UCHII regions (i.e., log(F6o/Fi2) ^ 
1.30 and log(Frj/Fi?) > 0.57). However, most of the 



Churchw ell et al.l (ll990) sources, and all of the lKurtz & Hofner 
d2005|) sourcesThave IRAS 100 11m > 1000 Jy and IRAS 
60 fim > 100 Jy and are therefore much m ore b iased towards 
brighter sources than the IPalla et al.l (119911) and ICodella et al 
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Table 3. Summary of number of sources observed, maser de- 
tections and detections rates by source classification. 



Source type 


# of Obs. 


# of Association 


Detection rate 


YSOs 


233 


60 


-25 


±4% 


HII/YSOs 


105 


34 


-32 


±5% 


Young/Old 


40 


7 


-18 


±7% 


HII regions 


125 


34 


-27 


±4% 



d 1995b surveys. Additio nally, the IChurchwell et alj (1990) and 
Kurtz & Hofnerl ( 12005b surveys are approximately 10 times 
more sensitive than the other two surveys. If a similar sen- 
sitivity to that of Pall a an d Codella were applie d to the 
Churchwell et all dl990h and iKurtz & Hofnerl d2005l) surveys, 
their detection rates would be reduced to 35% and 36%, re- 



spectively (see IKurtz & Hofnerll2005l for more details). Given 
the differences in sensitivity, source selection of these surveys, 
and statistical noise we find our results are in good agreement 

with all o f the surveys ment i oned. 

While lChurchwell et alj dl990h and lKurtz & Hofnerl d2005l) 
concentrated on the UCHII sta ge (identified from their ra- 



dio continuum em ission), the iPalla et al.l (119911) and the 



Codella et alj (119951) samples are likely to contain a mixture of 
nearby low-mass YSOs and more distant HII regions and mas- 
sive YSOs. Having a well selected sample of bona fide YSOs 
and UCHII regions associated with water masers we can now 
go a step further and look at how the association rates vary be- 
tween the different evolutionary stages. In Table [3] we present 
a breakdown of the number of RMS-H2O maser associations 
as a function of source classification (the errors are calculated 
assuming a binomial distribution). We can immediately see 
that the detection rates of H2O masers for all UCHII regions 
and YSOs are very similar. This would suggest that the con- 
ditions necessary for maser activity are equally likely in these 
two stages of the star formation process. The lower detection 
rate seen towards the Young/old class of sources is probably 
a reflection of the number of evolved stars within this group 
as discussed in Sect. 12.11 however, given the small number of 
statistics this may not be all that significant. 

It is somewhat surprising to find the detection rate for 
UCHII regions and YSOs is effectively the same, given that 
one might expect the environments to be very different. Water 
masers are widely thought to be associated with molecular out- 
flows and/or accretion, however, both of these phenomenon are 
expected to be disrupted and eventually halted altogether once 
the HII region begins to form. We compared the measured pa- 
rameters (e.g., maser velocity dispersions, peak intensity) for 
the UCHII regions and MYSO candidates to try and identify 
any differences that might give an insight into their nature, 
however, we found no significant differences - Kolmogorov- 
Smirnov (KS) tests of these parameters are consistent with the 
masers being drawn from the same overall population. 

The detection rate found tow ards YSOs and U CHII re- 
gions is double the rate reported bv lWang et al. I d2006l) (-12%) 
from observations of 140 compact cores found within a sam- 
ple of infrared dark clouds (IRDCs). These observations were 
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Fig. 5. Histograms of the distance from the Galactic centre of 
all YSOs and HII regions (outlined by solid line) and those 
associated with H2O masers (filled histogram). The bin size is 
1 kpc. 
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Fig. 6. Plot of the detection rate of H2O masers as a function of 
distance from the Galactic centre of all YSOs and HII regions. 
The bin size is 1 kpc. 



made with the VLA and had a sensitivity of -0.1 Jy beam -1 
channel -1 , which is a factor of a few times more sensitive than 
the observations presented here. The rate of detection of water 
masers is therefore significantly higher for star forming cores 
than for cores found in IRDCs, which are presumably at an ear- 
lier stage in their evolution, and many of which, are likely to be 
of relatively low luminosity. 



4.3. Galactic distribution 

Comparing the Galactic longitude and latitude distributions we 
find no significant differences between the RMS sources asso- 
ciated with masers and those without. However, we do find a 
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difference in the proximity to the Galactic centre between the 
two samples, with a much higher proportion of those associ- 
ated with H2O masers being located at smaller Galactocentric 
radii than those that are not associated. In Fig. [5] we present 
a histogram of the distribution of Galactocentric radii for the 
whole sample (outlined by the solid line) and the RMS-H2O 
maser associated sources (filled histogram). The distributions 
of the two populations are similar with both showing the same 
features. The two most interesting features are the two peaks 
located at ~6-7 kpc and 8-9 kpc. These distributions are very 
similar to the radial distribut ion of an IRAS selecte d sample 
of UCHII regions reported by lBronfman et al. I d2000h from CS 
observations. 

The 6-7 kpc peak correlates with a peak in the radial dis- 
tribution of a sample of southern infrared dark clouds (IRDCs) 
reported by Ijackson et al. ( 2008 ). This peak is at a larger ra- 
dial distance from the Galactic centre than found for a sam- 
ple of northern IRDCs which peaks at ~5 kpc. The differ- 
ence galactocentric distributions betwee n the Galactic first and 
fouth quadrants led ljackson et alJ (120081) to conclude that these 
features are more likely to be associated with the Scutum- 
Centaurus arm in a two arm G alactic model (e.g., Drimmell 
2000t iDrimmel & Spergelll2001l) t han part of a molecu lar ring 
of material located at ~5 kpc (see lJackson et al. 2006 and ref- 
erences therein). The coincidence of the peak in the distribu- 
tion of our sample of UCHII regions and M YSO candidates. 



the sa mple of IRAS selected UCHII regi ons (iBronfman et al 



2000l) and IRDCs djackson et alJ 120081) with the Scutum 



Centaurus arm is consistent with observations of external spiral 
galaxies where the formation of OB stars is seen to be exclu- 
sively associated with the spiral arms. 

The second peak located ~8-9 kpc is most likely due to a 
high number of sources found in the solar neighbourhood (^ 
0.5 kpc from the Sun). The distribution of the RMS sources 
associated with H2O masers shares some similarities with the 
distribution of the general RMS population, however, compar- 
ing it with the distribution of unassociated RMS sources reveals 
them to be significantly different. A KS test of these samples 
results in only a 6% probability that they are drawn from the 
same population. 

Looking in a little more detail we find that the detection 
rate varies as a function of distance from the Galactic cen- 
tre. This is illustrated in Fig. [6] which shows the detection rate 
to be approximately constant (~27%) between 4-6 kpc, after 
which it increases to ~37% between 6-7 kpc, before declining 
steadily to ~24% at 1 1 kpc. After 1 1 kpc the detection rate falls 
off to ~15%, however, the low number statistics for distances 
greater than this result in the spike between 12 an d 13 kpc. A 



long te rm study of water masers towards YSOs by lBrand et al 



(2003) found the detection rate, and maser stability, was sig- 
nificantly higher for more luminous young stars (> 3 x 10 4 
L Q ) than for lower luminosity YSOs. The increase in detec- 
tion rates for our sample of young massive stars located within 
the solar circle probably reflects the presence of a higher pro- 
portion of more luminous sources. Since the majority of our 
sources are located within the solar circle, and are thus subject 
to the distance ambiguity problem, reliable distances are not 



yet available for many sources. We are therefore not currently 
in a position to investigate this inference in more detail. 

In Fig. [7] we present a plot of the distr ibution of ^CO a s 
a function of Galactic longitude and Vlsr (IDame et al.ll200l[) . 
We have over-plotted the velocities of all detected masers; we 
distinguish between masers associated with an RMS sources 
and the unassociated masers by red and blue crosses respec- 
tively. The image presented in Fig. [7] illustrate the excellent 
correlation of RMS-H2O maser associations with the molecular 
material in the Galaxy as traced by the integrated 12 CO emis- 
sion. There is also a reasonable level of correlation between 
th e maser velocities and th e sp iral arm veloci ties as modelled 
by lTaylor & Cordesld 1993b and lCordesld2004 . 

In total 27 H2O masers were not associated with an RMS 
source (i.e., the velocity of their peak component is outside 
than 15 km s _1 criteria required for an association), this cor- 
responds to of all the detected masers. These unassoci- 
ated masers are the result of chance alignments along the line of 
sight. The high number of unassociated masers is due to the low 
resolution of these observations and the density of star form- 
ing regions located in the inner Galaxy. These masers fall into 
two catagories; those with velocities within the velocity range 



covered by the molecular line observations (see lUrquhart et al 



2007bl l2008al for details), and those at velocities outside the 
range normally associated with molecular gas in the Galaxy. 
We find weak CO emission at similar velocities for a few of the 
masers that have velocities within the Galactic range, however, 
for the majority no such emission has been detected. This is 
probably due to the difference in resolution of the CO observa- 
tions and the maser observations presented in this paper (which 
are ~34" and ~2' respectively), with the star forming regions 
associated with the maser lying outside the CO beam. 

The second category include masers at velocities outside 
those typically associated with the Galactic molecular gas. 
These can be clearly seen in Fig.|7]to have significantly differ- 
ent velocities (e.g., Z=10.4° and V L sr=72.7 km s" 1 and /=16.8° 
and Vlsr=70.4 km s~'). Since the conditions required to excite 
water maser emission are high densities (> 10 8 cm 4 ) and tem- 
peratures up to 400 K the lack of any large scale molecular gas 
is puzzling. One possibility is that these masers are associated 
with evolved stars, however, the nature of these masers remains 
uncertain and further observations will be required to ascertain 
their origin. 

5. Summary and conclusions 

We have presented the result of a set of water maser observa- 
tions towards a sample of 499 massive star forming regions. We 
detected a total of 163 masers above a Act sensitivity limit of 
~1 Jy, which makes this one of the most sensitive water masers 
surveys so far conducted. The majority of masers detected were 
previously unknown (-75%). We present plots for all detected 
H2O masers and tabulate their properties. 

By comparing the peak intensity velocities of the detected 
H2O masers with the velocities of the RMS sources we have 
identified 135 RMS-H2O maser associations. The overall de- 
tection rate towards our sample of young massive star forming 
regions is ~27%. This is in agreement with other water maser 




Galactic Longitude 



Fig. 7. Galactic longitude-velocity plot showing the velocities of all H2O masers detected as a function of Galactic longitude. 
Masers associated with an RMS source are s hown in red, whil e unassociated H2O masers are shown in blue. The distribution 
of molec ular material is shown in grey scale ( Dame et al. 200 lb for comparison. The location of the spiral arms taken from the 
model bv lTavlor & Cordesl (11993b and updated bv lCordes 



d2004 ) are over-plotted in yellow. 



surveys taking into account their various selection criteria and 
sensitivities. We find similar detection rates for MYSO can- 
didates and UCHII regions suggesting that maser activity is 
equally possible during these two stages of the star formation 
process. The detection rate as a function of distance from the 
Galactic centre is significantly enhanced within the solar cir- 
cle, peaking at ~37% ~ 6 kpc, possibly indicating the presence 
of a high proportion of more luminous YSOs and HII regions. 
Comparing our results with those reported by previous surveys 
towards massive star forming regions we find them to be con- 
sistent. 

These observations are a first step in our programme of 
follow-up observations designed to examine the global char- 
acteristics of this galaxy-wide sample of massive young stars. 
Further observations are planned to follow up the H2O masers 
detected from the observations presented in this paper the re- 
sults of which will be the focus of a subsequent paper. 
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